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Abstract 

A two-dimensional lattice system of non-interacting electrons in a homoge- 
neous magnetic field with half a flux quantum per plaquette and a random 
potential is considered. For the large scale behavior a supersymmetric theory 
with collective fields is constructed and studied within saddle point approxi- 
mation and fluctuations. The model is characterized by a broken supersym- 
metry indicating that only the fermion collective field becomes delocalized 
whereas the boson field is exponentially localized. Power counting for the 
fluctuation terms suggests that the interactions between delocalized fluctua- 
tions are irrelevant. Several quasi-scaling regimes, separated by large cross- 
over lengths, are found with effective exponents v for the localization length 
In the asymptotic regime there is v = 1/2 in agreement with an earlier 
calculation of Affleck and one by Ludwig et al. for finite density of states. 
The effective exponent, relevant for physical system, is v = 1 where the co- 
efficient of is growing with randomness. This is in agreement with recent 

high precision measurements on Si MOSFET and AlGaAs/GaAs samples. 
PACS numbers: 73.40.Hm, 71.55.Jv, 72.15.Rn, 73.20.Jc 
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I. INTRODUCTION 



The transition between quantum Hall plateaux in a two-dimensional electron gas is 
characterized by a divergent localization length £j with a critical exponent v and a non-zero 
longitudinal conductivity a xx . £j is finite and cr^ is zero inside the Hall plateaux whereas 
a xy is a constant. 

A direct measurement of the localization length exponent in a AlGaAs/GaAs sample by 
Koch et al. gave a value for v very close to 7/3. Recent high precision measurements on 
Si MOSFET |^],|3| and AlGaAs/GaAs samples [ffl, however, indicate that diverges with 
the electron density n like ~ b n (n c — n)^ 1 or with the magnetic field if like ~ bH{H c — H)~ x , 
where the quantum Hall transition (QHT) is at n = n c or H = H c , respectively. The 
exponent v « 1 appears to be almost independent of the material or the Hall plateaux. On 
the other hand, the coefficients b n , bn are sensitive to disorder: they increase with increasing 
disorder [0,0. This is a remarkable observation because in the scaling theory of Anderson 
localization || the coefficient is related to the mean free path. That means it would decrease 
with increasing disorder. The observation of v pa 1 is in sharp contrast to the experiment by 
Koch et al. The disagreement was explained in by insufficient sample size in the earlier 
experiment. 

The localization length scale near the QHT was also studied intensively in a number of 
numerical simulations using the network model ||f7| and the lowest Landau level approxi- 
mation PJlOll. These calculations agree on the result that the critical exponent is v ~ 7/3. 

Concerning the exponent v there is a calculation by Affleck |TT| based on the 
U(2n)/U(n) x U(n) nonlinear sigma model with topological term in the replica limit n — » 
||12|| . He obtains v — 1/2. The same value was found for Dirac fermions with random vector 
potential if the average density is finite [|1^ . 



A reason for the impressive agreement of the numerical calculations, on the one hand, 
and the disagreement between experiments, numerical and analytic calculations on the other 
hand could be the sensitivity of the QHT to the size of the system and to the type of 
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disorder. In particular, it may be related to the existence of a large characteristic length 
scale depending on disorder. The existence of such a typical scale is also indicated by the 
numerical results due to the fact that there is a cross-over from the pure network model 
[y = 1) to the random network model {y as 7/3). 

The purpose of this article is to investigate the role of disorder induced length scales in 
a tight-binding model with strong magnetic field. The work is based on an effective super- 
symmetric field theory for Dirac fermions with a random mass which enables us to study 
large scale properties. The main results are: 

1. Spontaneous breaking of the supersymmetry. I.e., the effective field theory for the QHT 
is not a nonlinear sigma model. 

2. The scaling behavior of the localization length depends on the characteristic scale 
exp(7r/g), where g is the strength of disorder: If <C exp(7r/g) the effective exponent 
is v — 1 whereas for ^> exp(7r/g) the exponent is v — 1/2. 

3. There is a universal value for the conductivity a xx = e 2 /h7i. 

The article is organized as follows: After the definition of the model in Sect. II an 
effective field theory is constructed for the averaged Green's functions (Sect. III). This field 
theory includes the description of the conductivity according to Kubo's formula. Then a 
collective field representation is introduced in order to cover symmetry breaking effects (Sect. 
III. A). The latter are discussed using a saddle point approximation for the collective field 
(Sect. III.B). Gaussian fluctuations around the saddle points and corrections to Gaussian 
fluctuations are studied in Sects. III.C and III.D, respectively. Finally, the localization 
length (Sect. IV) and the conductivity (Sect. V) are evaluated. 

II. THE MODEL 

A lattice model is considered in this article, stressing the universality in terms of the 
electron density and the magnetic field observed in the experiment, to study the asymptotic 
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behavior of the localization length near the QHT. Starting point is a microscopic model 
for non-interacting electrons on a regular lattice in a homogeneous magnetic field. Disorder 
enters only through a random potential on the lattice. This choice guarantees that the 
disorder does not affect the homogeneity of the magnetic field. The model is defined by the 
tight-binding Hamiltonian on a square lattice with magnetic flux <fi = Ba 2 , where a is the 
lattice constant and B the homogeneous external magnetic field. There is nearest neighbor 
hopping with rate t and next nearest neighbor hopping with rate t'. The Hamiltonian reads 
in Landau gauge 

H = — Y\te 2ivBay/(i,0 c{r)c\r + e x ) + tc{r)c\r + e y ) + t'e 2i7rBaiy± ^ /<t,0 c(r)c^r + e x ± e y ) 

r 

+h.c.) + V(r)c(r)c\r). (2.1) 

r 

e X)V are lattice unit vectors, and c' and c are fermion creation and annihilation operators, 
respectively. V(r) is a random potential representing disorder on the lattice. Without 
disorder, i.e., for V(r) = 0, this model was discussed extensively in the literature [p^IB|. A 
central result is the occurence of electron bands with a quantized Hall conductivity in each 
band. Gaps can be created in the model, for instance, by choosing a staggered potential 
V(r) = (— l) ri+r2 /i in the Hamiltonian ( j2.1| ) [I3|. By varying the staggered chemical potential 



one varies the concentration of electrons in the system. There are other methods to create 
gaps in a tigh-binding model. E.g., one could vary the magnetic field. This, however, would 
lead to a more complicated situation because the corresponding vector potential depends 
on space. In general, the relevant parameter for the quantum Hall transition is the filling 
factor n$>o/B. This is essentially determined by the ratio of the concentration of electrons 
n and the magnetic field B. Therefore, the variation of the concentration of electrons (i.e., 
the chemical potential) is equivalent with variation of the magnetic field in the quantum 
Hall system. 

In general, the creation of new bands ( "gap opening" ) can be described by Dirac fermions 
n] |TE| . Starting from the tight-binding Hamiltonian the Dirac fermions can be derived in 



a large scale approximation. The simplest case is that with half a flux quantum per lattice 



plaquette (0 = 0o/2) [p~3] , p79|| . (Such a strong magnetic flux is unrealistic in real crystals 
but typical for arrays of quantum dots in moderate magnetic fields [^0|.) For half a flux 
quantum per plaquette it is easy to derive the Dirac theory from a sublattice representation 
which takes into account the phase factor e lny ^ a of the tight-binding Hamiltonian and the 
staggered potential. The Fourier components of the non-random part H(k) read 
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with ( = it' / '4. All elements of the matrix are measured in units of the nearest neighbor 
hopping rate t. After expansion of k = (±7r, ±7r) + ap for small p vectors around the four 
nodes and a global orthogonal transformation H (k) — > OH(k)0 with 

(T -icr ' 
zer -cr 



(2.3) 



the Hamiltonian becomes 
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(2.4) 



The corresponding Green's function 

'Hu + iuo H 12 
H 2 i H 2 2 + iuj 



-1 



G = 
(H 11 + zuj)-' s 



(H u + iu)- x H l2 {H22 + ioj)- 1 
{H22 + tu)- l H 21 {H ll + iuo)- 1 [H22 + ioj)' 1 

decays asymptotically into two diagonal blocks 

(Hn+tLo)- 1 

(H22 + H0)- 1 



(2.5) 



(2.6) 



The lattice constant a is implicitly scaled out in t', //' and pj. ~ means asymptotics with 
respect to fi + 4i( ~ and p x ~ p y ~ 0. Thus the approximation breaks up the Hamiltonian 



( p.ip into two independent Dirac Hamiltonians H n / 22 = cr -p + cr 3 (/i=pt') with Pauli matrices 
(jj. The two Dirac theories describe particles with different masses fi =F t', respectively. The 
next nearest neighbor hopping term lifts the degeneracy of the two Dirac particles. Therefore, 
it plays an important role in this model and must be taken into account. Variation of the 
chemical potential implies a variation of the Dirac mass. In particular, Dirac fermions 
undergo a Hall transition if the mass vanishes [|l3|,[n[] . This is a consequence of the fact that 
the mass breaks the time-reversal symmetry: depending on the sign of m there is a clockwise 
or counterclockwise Hall current. If the light Dirac particle undergoes a Hall transition at 
fi — t' = its contribution to the Hall conductivity changes from a xy = — 1/2 for fi — t' < 
to a xy = 1/2 for fi — t' > 0. (The conductivity is in units of e 2 /h.) The heavy Dirac particle 
contributes a xy = 1/2 because its mass is positive. Thus the combined effect is a Hall step 
from a xy = to o xy = 1. This picture is particularly simple for $ = $0/2 but should also 
hold for other values of the flux as long as the low energy excitations are linear and can be 
described by Dirac fermions. 

For large scale properties like for the critical behavior near the Hall transition it is 
sufficient to consider only the light particle with // — t' 

(H' + iu)- l =[ =G{iu)=\ , (2.7) 

V*Vi-V 2 iuj-ix + t'j \G 2l G 22 J 

where V is the lattice gradient operator. Disorder, originally introduced in H by the random 

potential V, appears in Hu/ 22 as a diagonal matrix V with independent random elements 

V\, V 2 . The appearence of two random variables per site is a consequence of the sublattice 

representation required by the phase factor e my l a of the tight-binding Hamiltonian. The 

random matrix V is equivalent to a random mass Sfia^ and a random energy SEa^. For 

technical reasons the random energy term will be neglected in the following. 

It should be noticed that the random mass is marginally irrelevant on a perturbative 



level ||22j| . However, going beyond perturbation theory, it turns out that the random mass 



leads to spontaneous symmetry breaking which changes the properties significantly J|3 . 



This effect has not been included in previous studies of the localization properties near the 
Hall transition. It will be important for the considerations in this article. 

Dirac fermions can also be derived as the large scale approximation of the network model 
23 1 . It was recently pointed out by Ho and Chalker that the network implies a random 



Dirac mass (due to fluctuations in the tunneling rates), a random energy (due to fluctuations 
in the flux per plaquette) and a random vector potential (due to fluctuations in the phase of 
the hopping elements). I.e., in terms of the network model the random Dirac mass requires 
a fixed flux per plaquette and a fixed phase for the current between the vertices of the 
network. This is probably the simplest situation for the realization of a QHT. 

After averaging with respect to the random mass the localization length measured 



in lattice units a = ^J(f) /2B for electrons in a magnetic field B, is defined as the de- 
cay length of the function Cjf(r,u) = (\Gjf(r, 0; iu)\ 2 ). The relation \Gjf(r, 0; iu)\ 2 = 
Gjj/(r, 0; iu)Gjij(0, r; — iu) means that £j is given by the product of two Green's function at 
frequencies with opposite sign (retarded and advanced Green's functions). Due to the 2x2 
block structure of G there exists a relation between Green's functions at iu and Green's 
functions at — iu: 

Gjj(r,r'; -iu) = -Gy^r', r; iu), G jr (r,r'; -iu) = -G jf (r',r;iu) (j ^ j'). (2.8) 
This identity reflects the Lorentz-covariance of the Dirac theory. It implies 

\Gjj(r, r'; iu)\ 2 = —Gjj(r,r';iu)Gj>j'(r,r';iu) (2.9) 

and 

\Gjf(r, r'; iu)\ 2 = —Gjf(r, r'; iu)Gfj(r, r'\ iu). (2-10) 

This means that only the Green's functions with one frequency is required for the evaluation 
of localization properties in the relativistic model. The averaged quantity Cjj>(r, u) is trans- 
lational invariant. Therefore, it can be expressed by its Fourier components Cjj'(k,u). This 



can be used to calculate the localization length £j. (The following discussion holds for any 
choice of j, j'. Therefore, these labels are not written explicitly.) The correlation function 
C(r,u) for large r is proportional to r _a exp(— with some exponent a for which we 
assume that it is fixed for the model and does not depend on the parameters. This implies 
that 

"E r r 2 C(r,u) = E r r 2 -"exp(-r/^) 2 £ x x 2 ~" exp(-x) 
E r C(r,uj) ~~ £ r r- a exp(-r/&) " E x x~ a exp(-x) ' 

Dropping the constant term from the ratio of sums on the right hand side of ( 2.111) , the 

localization length can be defined in terms of the Fourier components as 



VlC(k,u) 



(2.12) 



\ C(k,u) 

The localization length is finite for w ^ but diverges in the regime of delocalized states 
with ijj — > 0. 

The localization length for massless Dirac fermions without disorder diverges like \io\~ x if 
uj = is approached. This behavior is probably unstable against arbitrarily weak random- 
ness, as it will be shown in this article. However, it has been shown in a previous paper 
that the localization length of the averaged correlation function has a lower bound which is 
the energy-energy correlation length of the 2D random bond Ising model. Since the latter 
diverges at the critical points, this implies a divergent localization length at the QHT. 

The calculation for the network model of Chalker and Coddington || , indicates that the 
critical exponent v = 1 of the pure model may change to v — 7/3 due to disorder. This 
was discussed as a possible appearence of a new random fixed point of the random model 
13fl . However, a new fixed point with such a behavior has not been found so far in terms of 



renormalization group calculations. 



III. FUNCTIONAL INTEGRAL REPRESENTATION 



It is convenient to introduce a functional integral representation for C(r, u), because this 
provides a basis to apply an approximation using a saddle point integration. The product 
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of Green's functions on the right hand side of ( |2.9|) and ( 2.1 0|) can formally be written as 



(iu + H);l r , f (iu + H T );,\ f)rk = J Xr'fXr^rk^r'k' exp(-S )V^V^V X Vx (3.1) 
with the quadratic form of the superfield (x r , ^r) 

_/Xr\ (iu + H \ fXr'\ 

S = -isign(u) £ ■ (3.2) 

r,r> V*J V iu + H T J . \y r r J 



X is a complex field and \l/ a Grassmann field, respectively. It is important to notice that H 
appears in the quadratic form for the complex field whereas H T appears for the Grassmann 
field. This difference will turn out to be crucial for the localization properties of the Dirac 
fermions. In particular, it will give all the delocalized states expected near the Hall transi- 
tion. In contrast, the quadratic form where H is used in the Grassmann sector instead of 
H T does not give these critical properties [ 26|, 21]] . 



Averaging with respect to disorder leads to 
((iu + H)^ rff (iu + H T );,l frk ) = J Xr'fXr^rk^r'k'(eM-So))V^V^VxVx. (3.3) 

A Gaussian distribution of the random Dirac mass is assumed in the following with mean 
m and variance g. Then the average can be performed exactly giving an additional quartic 
interaction term in So with coupling constant g. For weak disorder, i.e., small g, one could 
apply perturbation theory. Unfortunately, this does not lead to interesting results because 
it can not catch spontaneous symmetry breaking. In order to deal with the latter one must 
construct a representation which describes the field which is the conjugate to the symmetry 
breaking terms of the Dirac theory, the mass m and the frequency u. The appearence of the 
symmetry breaking terms in So dictates to choose XrXr and ^/ r ^/ r as the collective fields. 

A. Collective Field Representation 

In general, products the fields Xi ^ i n $o can be replaced by the collective fields as 
XrXr —* Qr, ^ r ^r —> P r , Xr^r @r and ^ r Xr ©r- (Some care is necessary to choose the 
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paths of integration correctly p7| .) One obtains for (|3.3| ) in collective field representation 
(for details see [ f27|j ) 

\ [e r>kj e r/>fk/ exp{-S')VPVQVeve {r^r') (3.4) 
9 

with the supersymmetric effective action 
S' = - ]T(Tr 2 Q r 2 + Tr 2 P 2 + 2Tr 2 e r 9 r ) 

9 r 

+ In det [{Hq + iuxjQ — iujctq + 2irPr) ] 

+ In det [1 - 4r6r(if T + iua + 2itPt)- 1 tQt(H q + iua - 2tQt)~ 1 ] . (3.5) 

Hq = ia ■ V + mo"3 is the average Dirac Hamiltonian and r the diagonal matrix (l,z). 
The introduction of the collective fields is important to discover the finite length scale e n ^ g , 



created by disorder, which is crucial for the properties of the random Dirac mass model [27 



B. Saddle Point Approximation 

A saddle point (SP) approximation is a crude approach for a two-dimensional system 
because it usually gives wrong results for low dimensional systems due to strong effects of 
fluctuations. In the model under consideration it will be used as a starting point to study 
also the fluctuations. An argument in favor of a SP approach is the fact that some features 
of the model can be described which are not available from pertubation theory. An example 
is the creation of states in the massive Dirac theory due to randomness [p3p8| . The hope 
is that the fluctuations around the SP are controlled by Gaussian fluctuations, at least if 
randomness is weak. This will be supported by the discussion presented below. Another 
argument for the SP approximation is its equivalence with the N — > oo-limit, where iV is a 



formal extension of the model to one for electrons with iV states per lattice site |f2~T |. 

The SP of the functional integral is given by the equations 5qS' = 5pS' = 0. The two 
SP equations are identical if one substitutes P = iQ: 

a 3 TQ r Tcr 3 = g(H + iua - 2tQt)~1. (3.6) 
10 



An ansatz for a uniform SP solution reads tQqt = —(irjao + m s <j^)/2. The SP equations 
imply a shift of the frequency uj — > rf = 77 + u; with 

rf -u = rj'gl (3.7) 

and a shift of the average mass m — > m' = m + m s with 

m s = -mgl/(l + gl), (3.8) 

7 ~ - / (r/ 2 + (m + m s ) 2 + k 2 )~ l kdk ~ — - ln[r/' 2 + (m + m s ) 2 ] = — In (3.9) 

7T JO 27T 7T 

with /i = m + m s + 277'. In the pure limit g — > the SP equations lead to 7/ = and 
m s = 0. For a given g > the SP depends on two parameters, m and cj. For u = and 
large |m| there is only a trivial solution of (|3.7|) with 77' — 77 — because gl < 1. As one 
varies |m| there is a critical point m c = 2/i c = 2e _7r / 9 where gl < 1 approaches gl = 1. As 
a consequence, the SP solution of (|3.7|) bifurcates from 77 = to 77 7^ at ji = /j c , and 77 = 



becomes unstable ||27|| . In the following only the region with 77 7^ is considered, where 



sign(rj) = sign(uj). This has a non-zero density of states 

p(m) » (l/27rp) limT/' = (l/2-Kg)rj = (l/A-Kg)J(m 2 c - m 2 )9(m 2 - m 2 ). (3.10) 

which describes a semicircular behavior. Of course, this must be normalized with an energy 
cut-off dependent constant. 

As an ansatz for the SP approximation with uj 7^ one can write 

r/M 2 = 7/(0, = 0) 2 + 6 = + 5j (3.11) 

where the last equation follows from the SP equation for uj = 0. This implies 

\p\ 2 = m 2 c /A + 5 (3.12) 

and for eq. (|3.7|) 

(g/2n)rf ln(l + 45/m 2 ) = w. (3.13) 
11 



The expansion of the logarithm for 5 <C m 2 c /A yields a cubic equation for 5 



,m 2 — m 2 _ x „ ,7im 2 



+ 5)5 2 « lu 



2 ^ , .2/ " '"c \2 



(3.14) 



Although this equation could be solved directly it is simpler to distinguish two different 

asymptotic regimes: 

(i) [m\ - m 2 )/4 < 5 < m 2 /4: 



and 

(ii) (m 2 -m 2 )/4 > 5: 



1 2<? ' 



(3.15) 



irm 2 1 
~ oo- 



9 \lm 2 — m 2 



(3.16) 



C. Gaussian Fluctuations 



In order to evaluate the localization length £ z the Gaussian fluctuations around the SP 
must be calculated (semi-classical approximation). Since Q, P and 6 are 2x2 matrices, 
the fluctuations can also be written as 4-component vector fields: q\ = 5Qu, q2 = (5Qu + 
SQ2i)/2, q% = —i(SQ\2 — 5Q2i)/2, q^ = 5Q22 with analogous definitions for p±, ...,p4 and 
for the Grassmann field ipi, ...,^4. The action of the Gaussian fluctuations reads in Fourier 
representation 

f 4 

S ' ~ J 2 [( l k)n,n'(qk,»q-k,v,> + Pk,nP-k,»>) + 2(l' k )^ip k ^ip_ k ^}d 2 k (3.17) 

with the fluctuation matrices l k and V k . For the large scale properties one needs only the 
asymptotic behavior for small wave vectors k. In particular, for a vanishing wave vector 
there is 

/l/#-2a/i* 2 2/5 

2/g-4a\fx\ 2 

2/g-4a\fx\ 2 

V 2(3 l/g-2afi 2 J 

12 
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(3.18) 
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fl/g-2afi* 2 

2(l/g-2a\fx\ 2 


\ o o 





2(3) 

2(l/g-2a\p\ 2 + 2(3) 










(3.19) 



l/g-2afi 2 J 



with 



a 



J (|/i| 2 + k 2 )- 2 d 2 k/Art 2 ~ \^\- 2 /Atx = l/nm 2 c 



(3.20) 



and 



(3 = J k 2 (\/j\ 2 + k 2 )~ 2 d 2 k/4ir 2 = 1/2 - l/4?r ~ - In |/x| 2 /4tt ~ l/2g. 



(3.21) 



These quantities become quite large for small indicating a short range behavior of the 
related modes. However, the second diagonal element of I' vanishes for vanishing to because 
of 1/g — 2a\fi\ 2 — 2(3 = w/gr]'. This is a direct consequence of the SP condition (|3.7| ) and 
indicates a critical mode for all \m\ < m c . Moreover, it implies a divergent behavior 
of the localization length if u — > 0. The corresponding correlation function of the critical 
mode can be calculated in the large scale limit by expanding Vk in powers of k as 



C(k,u) = g- 2 {2u/gr]' + Dk 2 ) 



2\-l 



(3.22) 



where 



D = Aa 



l + a{ 



+ 



+2 



(3.23) 



l/g-2afi 2 l/g-2afi* 2 ' 

grj'D/2 is like a diffusion coefficient. It is real and it never becomes zero. The critical behav- 
ior describes a phase with a divergent sum J2 r ,j(\Gjj(r,r')\ 2 ); i.e., the correlation function 
decays non-exponentially. It implies that Cik = 0,uj) = v( /2guj = n (p(uo)) / u . This holds 
not only on the SP level but in general due to the identity J2 r ',j,j' \Gjj'{j,r')\ 2 = ixpfuj. 

C in ( |3.22| ) is an approximation based on Dk 2 1. Since D diverges like 1/nm 2 . it is 
not possible to use this expression for the pure limit. This also reflects the non-perturbative 
character of the SP approximation. 
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Apart from the critical (delocalized) fermion mode there is a boson mode which becomes 
critical at m — ±m c . This is due to a vanishing eigenvalue of Iq because rj = at m = ±m c . 
Thus there are delocalized states for \m\ < m c due to massless fermions whereas a combi- 
nation of critical fermions and critical bosons controls the QHT. The band of delocalized 
fermions simplifies the study of transport properties of the model away from the critical 
points m = ±m c . 

An analogous calculation for the Gaussian fluctuations can be performed for Anderson 
localization (i.e., for <fi = in eq. Q2.1|) ). In that case the Green's function of eq. ( |2.7|) must 
be replaced by 

'iuo + M + V 2 

iu - M - V 2 

Now the matrix H = M + V 2 is symmetric in contrast to H = ia ■ V + (// — t')a^. As a 
consequence the corresponding fluctuation matrix is degenerated for the fermion and the 
boson sector (i.e., I' k = Ik) with 

fl/g-a 



G{iu) 



(3.24) 









2(1/9 -P) 

2(1/(7-/3) 








l/g-a*) 



(3.25) 



where 



a 



(k 2 -E + i7]) 2 
{27i) d J {{k 2 - E) 2 + r) 2 ) 2 



<fk 



^3.26) 



and 



-d d k. 



(3.27) 



(2n) d J (k 2 -E) 2 + 7] 2 ' 

There are two vanishing eigenvalues of Io for vanishing uj in the Grassmann as well as in 
the complex sector, each of them is 1/ g — (3 = uj/grj'. Therefore, the large scale behavior of 
this model is very different, and the critical properties are described by a nonlinear sigma 



model including fermion as well as boson degrees of freedom [29]. This is a consequence of 
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the fact that = preserves the supersymmetry whereas the model with = 0o/2 breaks 
the supersymmetry implying I ^ I'. The latter implies that there is only one massless 
(delocalized) fermion field, and all other fields are massive (localized). 

D. Corrections to Gaussian Fluctuations 

Gaussian fluctuations are usually not sufficient to describe the properties of a critical 
system, especially at low dimensionality, because interaction between the fluctuations are 
a relevant perturbation. For Anderson localization the interaction is marginal in d = 2 as 
one finds from power counting of the scaling behavior. This method can also be applied to 
the effective field theory of random Dirac fermions. It provides a first check for the effect of 
the interaction among the fluctuations. In the following we will see that the perturbation 
term for random Dirac fermions has dimensionality —2. Therefore, the interaction of the 
fluctuations is irrelevant. (It seems that the dimensionality is reduced by 2 in comparison to 
Anderson localization. This is similar to the dimensional reduction by 2 in supersymmetric 
theories, applied to the average density of states for a particle in a strong magnetic field and 
a random potential |30|j31[| .) 

Away from the critical points \m\ = m c it is sufficient to study the Grassmann fluctu- 
ations because the complex fluctuations are massive. Their action reads for terms up to 
fourth order in the fluctuations 

S"« 2 f(2u/g7 1 ' + Dk 2 )i) k ^. K2 d 2 k - 8 £ ^[e^G^e^G^^G^Qr^^], 

ri,...,r4 

(3.28) 

where G = tGqt. is the only critical mode of the collective Grassmann field O r . Thus 
the interaction term can also be written as 

-8 ^2^i(5 rir2 Cr 1 (5^ 2r3 Cr 1 (5 r 3 r4 (7i(5^ r J^ rij 2V'r 2 ,2^r3 ) 2V , r4 ) 2- (3.29) 

n,...,r 4 

The trace term can be evaluated and yields after some straightforward calculation together 
with the approximation G r y ~ \fi\~ 2 [(ir]' <jq — mas)S r y + iVi^y&i +i^2;r,r' a 2] the following 
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expression 

S|/x| [4 , r 1 dr 1 ,r2'4 , r2dr2,r 3 ' i Pr 3 d r3ir4 1p r4 d r4jri + ^ri )r2 ^r 2 ^r 2 ,r 3 ^3 ^ r3j r 4 ^4*^4, rj (3.30) 

r 1 ,...,r 4 

where the index 2 of the Grassmann field has been dropped and d r ^> = (Vi- r , r ' + «V2 ; r/)- 
This result reflects the fact that terms with an odd number of V-operators cancels in ( |3.29|) 
and terms quadratic in V cancel each other because of the anticommutation rule of the 
Grassmann field: ip ri ip ri 4>r 3 ipr 3 + 4'r 1 i , r 3 4'r 3 i , r 1 = 0- Higher order terms with at most sec- 
ond order gradients disappear individually because they contain Grassmann fields at the 
same site. Simple power counting indicates that this interaction term has dimensionality 
—2. Therefore, the interaction is irrelevant in comparison with Gaussian fluctuations con- 
sidered in the previous section, and it scales quickly to zero under renormalization group 
transformations. 

The situation is different if we approach m = ±m c because the complex field also becomes 
critical. As a consequence, the corrections to Gaussian fluctuations are marginal rather than 
irrelevant then. It is possible that the localization length becomes finite at m = ±m c due 
to renormalization effects. This phenomenon near ±m c requires a separate treatment which 
will not be considered in this article. 

The irrelevance of the interaction terms is very special for the model under consideration. 
In similar two-dimensional models, like the Gross-Neveu model or the tight binding model 
without magnetic field (orthogonal nonlinear sigma model) or with weak magnetic field 
(unitary nonlinear sigma model), the interaction is always marginal (i.e., dimensionality of 
the interaction term is zero). 



IV. LOCALIZATION LENGTH 

According to the discussion of the Gaussian fluctuations in Sect. IIIC the leading large 
scale behavior is given by the correlation function of ( |3.22| ) . From the latter the localization 
length can be calculated, using eq. ( |2.12| ), as £j = (2Dgr]' /cu) 1 ^ 2 . Together with the SP 
results ( |3.11|) and (|3.23|) reads 
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&~^(^+*) 1/4 «/4 + *r^. (4.D 

The localization length diverges like uj~ 1 / 2 because D and rf remain non-zero for \m\ < 
m c according to the results of the SP approximation. To compare £j with numerical or 
experimental results it is important that D can be large (~ l/^l -2 ) arid 7/ can be small. 
Therefore, several quasi-scaling regimes exist as indicated by the graph in Fig. 1 which is 
simply the plot of (|4.1|) together with ( |3.13|) . Only in the asymptotic regime u ~ 0, i.e., for 
[m 2 — m 2 )/4 ^> 5, the localization length diverges like u^ 1 ^ 2 

^-^r': m2 ) 1/4 -^- i/2 - (4.2) 



7r V 4 ' m c 

Surprisingly, the exponent v = 1/2 agrees with that of a completely different approach to 



the QHT by Affleck [IT]. Moreover, u = 1/2 was also found by Ludwig et al. for Dirac 
fermions with a random vector potential at that special point on their fixed point line where 
the average density of states is finite 



A quasi-scaling regime occurs for (m 2 — m 2 )/4 <C 5 <C m 2 /4 where we have 

5 „ ,/^(m c 2 /4 + fl-Va^ )i/B w -i/3 „ ,&1(^)^-V3 (4.3) 

V 7T 2(7 V T 771 c 2(7 

i.e., the effective exponent is z/ = 1/3. 

Going back to the SP equation ( [3.131 ) one could also assume m 2 /4 <C 5. It gives 

-5 1/2 (g/2n)\n(m 2 c ) = 5 1/2 (4.4) 

This result implies for the localization length 

£ „ &-i/4 w -va „ ( 4 . 5) 

V 7T V 71 

i.e., the effective exponent is v = 1. The coefficient in (|4.3|) decays for small (7 like g 1 / 3 e 2n ^ 3g 
whereas the corresponding coefficient in fl4.5|) grows like g 1 / 2 with the strength of disorder. 
The coefficients of the power law for v — 1/3 (ci), for v — 1/2 (02) and for u = 1 (03) are 
plotted in Fig. 2. z/ = 1 and the behavior of the corresponding coefficient C3 agree with the 
observation in experiments |fj [| . 
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V. CONDUCTIVITY 



The longitudinal conductivity can be calculated via Kubo's formula 

e 2 

<7xx{u) = -^J 2 r2 (Gjj'( r , 0; iu)Gjfj(0, r; -iu)). (5.1) 

The correlation function is again the expression we have considered in the effective field 
theory. Since for small us only the large scale part of the correlation contributes significantly, 
we can use C(k,u) of ( |3.22| ) to write 

a xx (cu) = - e luj 2 V 2 k C(k^)\ k=Q = ^u 2 g- 2 D{u)/grf)- 2 = jD V ' 2 . (5.2) 

For weak disorder we use the approximation D ~ 4a. Furthermore, for u ~ we get 
a = l/47r |/z| 2 ~ l/7rm 2 . Therefore, the conductivity reads in the dc-limit 

2 2 2 

a xx (cu = 0) ~ e m c -m Q( ^ 2 _ ^ ^ 



nh m 



For m = the dc-conductivity is independent of disorder 

a c xx (u = 0)~e 2 /nh. (5.4) 

This result is in agreement with a calculation for Dirac fermions with a random vector 
potential [pL3 1 . Thus there is a universal conductivity in d — 2 at the center of the band of 



extended states — m c < m < m c . 

The Hall conductivity can be calculated using a simpler field theory which does not break 
the supersymmetry ||21||. In units of e 2 /h it was found 



o~xy ~ 1/2 + sign(m) 1/2 — (l/7r)arctan(y m 2 c /m? — l)6(m c — m). (5.5) 

For m = the Hall conductivity is always e 2 /2h whereas for m ^ it depends on disorder 
via m c . The resistivity, calculated from the averaged conductivities (which is an additional 
approximation because one should actually evaluate the averaged resistivity) yields 

PL = <J(° C J + < y 2 ) « 0.9/i/e 2 . (5.6) 



This value agrees with an error of about ±10% with various experimental results |32] 



VI. DISCUSSION OF THE RESULTS 



For a given strength of randomness g there are three different regimes of the behavior of 
the localization length, depending on the value of 5. 5 is directly related to the frequency uo 
according to the equations ( |3.15|) , ( |3.16|) and ( |4.4| ), respectively. The corresponding power 
laws ~ c j UJ ~ Uj appear with effective exponents v\ = 1/3 and u 3 = 1 and with the asymp- 
totic exponent v-i = 1/2. The exponents and the corresponding coefficients Cj are shown 
in Table I, and the g dependence of the coefficients is plotted in Fig. 2. Ludwig et al. 
evaluated the exponent v for two-dimensional Dirac fermions with random vector potential 
with variance g a . In that case the exponent of the asymptotic localization length depends 
smoothly on randomness as v = 1/(1 + g a /7i) [13|. It should be noticed that the inequality 
v > 1 of Chayes et al. |33| for two-dimensional random systems does not apply to the prob- 



lem under consideration (cf also the discussion in Ref. |]13|). It is remarkable though that 
there is agreement of the asymptotic result v = 1/2 between the present calculation for the 
random Dirac mass, the random vector potential with finite average density of states (where 
g a /^ = 1) |~3| and the nonlinear sigma model with topological term ||11|| . 



The result v = 1/3 is not reliable because it appears close to the critical points |m| = m c , 
where the bosonic degrees of freedom become critical. Moreover, the asymptotic regime 
with v = 1/2 is not realistic: The typical width of the fluctuations ^fg in a sample is about 
10% of the hopping rate. I.e., m c = 2e~ n ^ g is immeasurably small. Therefore, only the 
regime 5 3> m^/4 would be accessible, and only the effective exponent v = 1, together with 



the coefficient y2g/n, is of practical relevance. Consequently, only the coefficient of the 
localization length is affected by randomness. This is in agreement with the observations 
in Refs. |§-§]]. However, the divergency of the localization length is controversal among 
the various experiments. In this context it would be interesting if the cross-over length, 
evaluated in the article as ~ exp(jr/g), can be observed experimentally. 
The density of states is non-zero near the QHT. This result is non-perturbative because 
p oc exp(— n/g), and it reflects spontaneous symmetry breaking [E3f. It agrees with a Monte 
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Carlo simulation for the network model by Lee and Wang [341 , with an exact diagonalization 



of a finite system |j3"5| , with a rigorous estimation |23| and with an exact calculation for 
Lorentzian disorder [[28[ . 

The longitudinal conductivity a xx (u>) is non-zero between m = —m c and m = m c , the 
transition region between two Hall plateaux. In an experiment it may not be possible to 
resolve the width of this band of delocalized states since m c is too small. Therefore, the 
width of the transition between the Hall plateaux will always be dominated by thermal 
broadening. Thus a power law for the width A at temperature T like A ~ T K [^] is a 
realistic ansatz. The conductivity at the QHT should agree with the value of a xx at m = 0, 
i.e., a xx = e 2 /fin. Converting this value into the corresponding value of p c xx gives pa 0.9h/e 2 
which is in good agreement with experiments. The universal value a xx = e 2 /hir agrees with 
the value found for random vector potential |13| but not with the result found from the 
numerical simulation of the network model, where a xx = e 2 /2h was found 0. This value 
was also obtained for the lowest Landau level projection ||. 

It seems that our results for the 2D Dirac fermions with random mass are in good 
agreement with the corresponding results for 2D Dirac fermions with random vector potential 
and with recent experiments. But there is disagreement with the results of the numerical 
simulation of the network model. This may be related to different types of randomness (e.g., 
strong randomness in the magnetic field) or to different geometries. It is also possible that 
there is a strong renormalization of the localization length near m = ±m c due to extended 
boson fields leading to an exponent v pa 7/3. This, however, would raise the question about 
the origin of the experimental value v pa 1. Moreover, the disagreement of the values of a xx 
for the network model on the one side and for Dirac fermions with m = on the other side 
cannot be explained by renormalization effects. 
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TABLES 



TABLE I. Exponents for the localization length and the corresponding coefficients. 



regime 


5 <C (m 2 c - m 2 )/4 


(ml - m 2 )/A < 5 < 


5 > m^/4 


V 


1/2 


1/3 


1 


coefficient 


C2 


Cl 


C3 
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FIGURES 

FIG. 1. Scaling of the localization length as a function of the frequency for disorder strength 
g = 0.2 and /z - t' = 0. 

FIG. 2. Coefficients cj of the power laws of for /x — t' = 0. 
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Fig. 1: 
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